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Abstract 


Two-particle correlations at small relative momenta give insight into 
the size of the emitting source. Of particular interest is the test of hydro- 
dynamic models which predict a universal, apparent decrease of the extent 
of the system with increasing transverse mass mr as a consequence of the 
strong radial flow in heavy-ion collisions at LHC energies. This contribu¬ 
tion presents a study of correlations for protons and A particles in Pb-Pb 
collisions at -^snn = 2.76 TeV measured with ALICE. The investigated 
particle species expand the experimental reach in m t due to their high 
rest mass. Residual impurities in the samples from misidentification and 
contributions from feed-down are corrected using data-driven techniques. 
Correlation functions are obtained in several centrality classes and i»t 
intervals at large itit and show the expected decrease in volume of the 
strongly interacting fireball for more peripheral collisions. We observe 
the decrease of the source size with mx, predicted by hydrodynamical 
calculations, out to (mx) = 2.18 GeV/c 2 . 
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1 Introduction 


Femtoscopy aims at measuring the size of the particle emitting source. Experi¬ 
mentally, this is achieved by studying two-particle correlations at small relative 
momenta. If the size of the system is large, the emissions of the two parti¬ 
cles will likely be substantially separated in space. Therefore, no interaction or 
symmetrization will take place and the momenta will be uncorrelated. If on 
the other hand - the extent of the source is small, the vicinity of the emanation 
points will trigger the species-specific interplay with a characteristic dependence 
of the correlation on the momentum difference. Given that the two-particle in¬ 
teraction is known, a source size can be inferred from the correlation function. 

No Coulomb or quantum-statistical effects take place in the pA system and 
the strong interaction parameters are sufficiently well known from, e. g., bubble 
chamber experiments. A typical reaction for studying the pA final-state inter¬ 
action is the production of A baryons by shooting a beam of negatively charged 
kaons on a hydrogen target as done with the Saclay bubble chamber 11. The A 
hyperon - created via K _ + p —>• A + 7t° - can subsequently scatter elastically 
off another proton. Measuring the cross-section of the elastic process differen¬ 
tially vs. the excess energy directly quantifies the attractive final-state strong 
interaction, which turns out to be comparable to the nucleon-nucleon one. 

Based on the model by Lednicky and Lyuboshits [2], the sensitivity of the 
pA correlation function to the volume of the hot and dense medium created 
in heavy-ion collisions was first explored in ;3 . The correlation was found to 
be affected in height and shape by the size of the source. Furthermore, the 
interdependence of the two momenta was shown to keep its susceptibility to a 
change in source size for radii larger than 4 fm - an advantage over the Coulomb- 
depleted proton-proton correlations. 

A particularly interesting subject is the investigation of the reaction dynam¬ 
ics of the strongly interacting matter generated in Pb-Pb collisions. The large 
pressure gradients give rise to an expansion of the medium. The resulting col¬ 
lective velocity competes with the thermal velocity ^/T/mx; consequently the 
tot dependence of the radii probes the dynamics of the source. At low tot, no 
strong collective motion is present and the particles will be correlated over the 
full extent of the source. At high tot in contrast, the radial flow introduces a 
positive correlation between the emission point and the particle’s momentum. 
Looking only at pairs of particles with a small momentum difference, the cou¬ 
ple’s constituents will originate from the same region, which is smaller than the 
geometrical size of the fireball, leading to the apparent shrinking of the particle 
source with tot- To probe the high tot regime, it is experimentally beneficial to 
investigate particles with a high rest mass. So far, the pA system is the heaviest 
studied system [4j|6]. 
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2 Particle Selection 


The analyzed Pb-Pb data was taken with ALICE 7 . The TO and VO scintil¬ 
lator arrays provided an event-trigger signal and an estimate for the collision 
centrality. The main detectors for particle tracking in ALICE are the Inner 
Tracking System, Time Projection Chamber (TPC), and Transition Radiation 
Detector. Particle identification for this analysis was performed by the TPC 
and the Time Of Flight (TOF) detector. The vanishing ^ at the LHC entails 
that matter and anti-matter are produced in equal abundance. While the corre¬ 
lation functions for pairs of pA and pA were obtained separately, the selections 
for particles discussed in the following hold true for their charge conjugates as 
well. 

The proton selection depended on the reconstructed momentum and used the 
TPC and/or TOF. The resolution of the truncated specific energy loss measure¬ 
ment {dE/dx) can be parametrized with a Gaussian of width a. The TPC pro¬ 
vided a separation of the most probable d E/dx measurement for protons to the 
one of any other species of more than four a up to p = 0.75 GeV/c. Above this 
kinematic restriction, the expected dE/dx for protons and the ultra-relativistic 
electrons becomes too similar, which impeded an unambiguous particle identi¬ 
fication by the TPC alone. For higher momenta up to p = 5.0 GeV/c, TOF 
allowed for an unambiguous proton identification |8j. With the chosen criteria a 
proton purity - defined as the number of protons over the number of all selected 
particles including misidentified tracks - above 99% was achieved over the full 
dynamic span. 

The A selection was based upon a V° topology finder, which reconstructs 
the charged decay A — > p + n~ ; the liyperon is identified via the invariant mass. 
Fig. a (left) shows a fit to the invariant mass distribution of all A candidates 
employing a Monte-Carlo template for the signal shape in an exemplary phase- 
space bin with a high yield, namely 0.5 < \y\ < 0.6 and 1.0 < pr (GeV/c) < 1.5, 
for the 10% most central events. The parametrization allows to determine the 
purity of the sample - defined as the ratio of picked A particles to all taken 
V° vertices - which amounts in this case to pur A = 91% within the window 
in invariant mass of ±4 MeV around the PDG value. Determining the purity 
enables one to correct the correlation function for the uncorrelated background 
in a following step. Fig. [l] (right) shows the evolution of the A purity as a 
function of rapidity and transverse momentum for the 10% most central events. 
Protons knocked-out from an interaction with the detector material resemble 
displaced tracks from weak particle decay vertices. A contribution of protons 
from material manifests itself in a degradation of the purity for pr < 0.5 GeV/c; 
the effect is absent for the case of A. For higher transverse momenta, the purity 
is 81% for \y\ > 1.0, within \rj\ < 0.9 it is better than 90%. 
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Figure 1: Left: Invariant mass spectrum of A candidates in an exemplary phase- 
space bin with a fit using a Monte-Carlo template for the signal shape. Right: 
A purity as a function of rapidity and transverse momentum for the 10% most 
central events. 


3 Feed-Down Determination 


The pA strong interaction is limited in its range to a few fm. Hence, any prod¬ 
uct from electro-magnetic or weak decays will not contribute to the excess seen 
at small relative momenta. The contamination in the proton sample mostly 
stems from the decay A —> Pdec7t~- Since no significant AA correlation was 
seen by the STAR Collaboration |[9], the decay proton will also not carry any 
residual correlation 10 ■ 12 from the mother particle^ The determination of the 
amount of this uncorrelated feed-down allows for a correction of the pA corre¬ 
lation function with the feed-down fraction / p = non^Mimjuy r^ e outstanding 
performance of ALICE allows to determine this feed-down fraction in the pro¬ 
ton sample directly from the data via the distance of closest approach (DCA) 
of the track extrapolation to the primary vertex. The two-dimensional (trans¬ 
verse and longitudinal) DCA distribution was obtained differentially in rapidity, 
transverse momentum, and centrality for the data, as well as for templates from 
Monte-Carlo simulations for primary protontj^J protons from weak decays, and 
protons from an interaction with the detector material. The distinct shapes of 
the templates - almost flat in DCA for the material contribution, wide for the 
weak decays, and peaked for the primaries - allow to disentangle the different 
origins of the charged nucleons by fitting the templates to the data, as pictured 
in Fig. [2] (left). Selecting only tracks with a DCA smaller than 1 mm in the 
transverse and smaller than 1.5 mm in the longitudinal direction enhances the 
primaries in the sample. In the exemplary phase-space bin 0.0 < y < 0.25 and 
1.0 < pt (GeV/c) < 1.5, the feed-down fraction / p totals to 15%. This arises 
from the amount of particles which are from weak decays and from material of 
15% and less than 1%, respectively. 


1 Recent data from ALICE |13| and STAR Jl4] suggest a slight AA anti-correlation. The 
momentum released in the A decay will wash this correlation out; the small fraction of p,] e(: A 
pairs in the pA sample will make it a tiny, likely negligible, correction. 

2 According to the common ALICE definition, primary protons include decay products, 
except products from weak decays of strange hadrons. 
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Figure 2: Left: Distribution of the transverse distance of closest of approach of 
protons to the primary vertex in data with a fit utilizing Monte-Carlo templates 
in an exemplary phase-space bin. Right: A pair purity for pairs with wt > 
1.9 GeV/c for the 10% most central events with a constant fit in the region 
0.0 < q (GeV/c) < 0.15. 


Also the A sample is contaminated by feed-down. The contributing, dom¬ 
inating weak decays are 2° —» Arc 0 (BR 99.5%), 2~ —► Arc - (BR 99.9%), 
—> AK _ (BR 67.8%), and Cl~ —> 2 (—> An)n (BR 32.2%). The reconstruc¬ 
tion efficiencies with the chosen A selection criteria in this analysis for all Y = 
2°, 2°, 2“, 2 + , n _ , and fl + were obtained three-dimensionally in p\ —>• (p^, y A ) 
for each centrality class of [l5] from a Monte-Carlo simulation. Using the p t- 
differential spectra of S _ , S + , f l~, and f l + measured by ALICE 
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for several 

centrality classes and assuming isospin symmetry for the unmeasured 2° and 
2°, enabled us to determine the fraction of A hyperons from weak decays. 

The cr = 22 pm of the electro-magnetic decay E° —>■ Ay (BR 100%) is 
much larger than the range of the strong interaction, but also too small to be 
resolved experimentally. In 16 it was found that the final-state interaction in 
the pZ° channel is much smaller than in the pA system. Thus, a significant 
fraction of A come from L° and are uncorrelated with primary protons. We use 
the fraction of A from L° properly taking into account all strongly decaying 
resonances - from the thermal model of 17 , while systematically considering 
the study in 18 with an additional variation in the freeze-out temperature, 
and the results of 19 with the value at the kinetic and chemical freeze-out. All 


models give that about 30% of the A originate from a electromagnetic decay. 
We unite the fraction from weak and electromagnetic hyperon decays in the 
feed-down fraction /a- 


5 


































4 Corrections and Results 


The raw pA correlation as a function of the momentum difference is obtained 
as the ratio of pairs from real events over those reconstructed from mixed 
events for three centrality divisions and up to four tot classes for each cen¬ 
trality class. Dealing with non-identical particles, we use the generalized mo¬ 
mentum difference introduced in 
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by R. Lednicky q = \q — P(qP)/P 2 \,q = 
Pi — P 2 , P — Pi + P 2 , where p± and P 2 are the momenta of the particles. In the 
following, we omit the tilde. The overall pair purity factorizes: pur(g, ?7 Jt) = 
Pur A (< 7 , m T ) • (1 - / p (<7,tot)) ■ (1 - /a(< 7, m T ))- The A pair purity pur A (g,m T ) 
is shown exemplary in Fig. [2] (right) for the 10% most central events and the 
highest ttit class, tot > 1.9 GeV/c 2 . Its value of 95.1% is constant over the 
region of interest q < 0.15 (GeV/c) and beyond. It is obtained by looking up 
the single-particle A purity in the ^differential histograms for the given 

centrality; the same holds for / p and / A accordingly. All impurities constitute 
- as discussed - an uncorrelated background. Hence, the corrected correlation 
function can be attained by scaling the raw correlation function with the inverse 
pair purity: 


c 2 corr (g,m T ) = 


1 


pur {q,m T ) 


(C™ w (q, m T ) — 1)) + 1- 


(1) 


The effect of the finite momentum resolution was studied with a Monte- 
Carlo simulation. For small relative momenta, i. e. (q rec ± q gen )/2 < 0.1 GeV/c, 
where q rec is the reconstructed and q gen is the generated relative momentum, the 
deviation, quantified as ( q rec — q S en)/V 2, can be parametrized with a Gaussian 
with a mean of —0.32 ±0.02 and —0.38 ±0.08 MeV/c and a width of 7.26 ±0.02 
and 7.13 ± 0.07 MeV/c for the 0-10 and 30-50% most central events. In an 
mT-differential study, the width turned out to slightly increase by less than 
2 MeV/c for higher tot • 

Fig. [3] shows a set of exemplary correlation functions corrected via Eq. [I] 
i. e. remedied for the uncorrelated background coming from weak and electro¬ 
magnetic decays. The correlation functions for pairs of particles and pairs of 


anti-particles were merged following the recipe of the Particle Data Group 21 


The systematic errors include a variation of the correction for the electro¬ 
magnetic feed-down, a change in the hyperon input spectra for the correction 
of the A from weak decays, an altered momentum resolution, normalization, 
invariant mass selection of the A, changed DCA cuts on the proton, a varied 
two-track resolution cut, and dominantly the uncertainty on the pA interaction 
when fitting the data. The top panel shows the centrality dependence, i. e. the 
0-10% most central events in red and the 30-50% most central events in blue; 
both samples have a mean transverse mass of 1.5 GeV/c 2 . One clearly sees the 
expected effect of an increased width and height of the excess at small relative 
momenta, which translates into a smaller source size, for the more peripheral 
collisions. The bottom panel shows the dependence on tot with the green points 
representing 1.0 < tot (GeV/c 2 ) < 1.4 resulting in (tot) = 1.27 GeV/c 2 and 
the orange symbols depicting tot > 1.9 GeV/c 2 yielding (tot) = 2.18 GeV/c 2 


6 





for the 0-10% most central events. Also here, a clear ordering is apparent with 
the higher to t giving evidence for a smaller source than at lower rax, matching 
the expectation within the hydrodynamic picture outlined in the introduction. 
Note that the correlation function with (tot) = 2.18 GeV/c 2 exceeds the (tot) 
of any previous measurement from 6| at the SPS, 5] at RHIC 
LHC. 
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5 Summary 

The supreme performance 
of ALICE makes it possible 
to collect very pure samples 
of protons and A with rich 
statistics. We obtained pA 
correlation functions multi- 
differentially in centrality 
and transverse mass. They 
were corrected for misiden- 
tification and contamina¬ 
tion from weak and electro¬ 
magnetic decays, employ¬ 
ing data-driven methods 
to quantify the impurities. 

The correlation functions 
presented here represent the 
largest tot reach of any 
femtoscopic measurement, 
with result being shown for 
to t > 1.9 GeV/c 2 . The 
conveyed centrality depen¬ 
dence of the pA correlations 
exhibits the expected be¬ 
havior of a smaller source 
for more peripheral colli¬ 
sions. The communicated 
TOx-differential correlations, 
spanning a range in (tot) of 

more than 0.9 GeV/c 2 , display the decrease in source size with tot, qualitatively 
agreeing with hydrodynamic predictions, out to highest (tot) = 2.18 GeV/c 2 . 
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Figure 3: Exemplary pA and pA correlation func¬ 
tion corrected for weak and electromagnetic de¬ 
cays. Centrality dependence at the top, Tot de¬ 
pendence at the bottom. 
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